Multi-band quantum ratchets 
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We investigate directed motion in non-adiabatically rocked ratchet systems sustaining few bands 
below the barrier. Upon restricting the dynamics to the lowest M bands, the total system-plus-bath 
Hamiltonian is mapped onto a discrete tight-binding model containing all the information both on 
the intra- and inter-well tunneling motion. A closed form for the current in the incoherent tunneling 
regime is obtained. In effective single-band ratchets, no current rectification occurs. We apply our 
theory to describe rectification effects in vortex quantum ratchets devices. Current reversals upon 
variation of the ac-field amplitude or frequency are predicted. 
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A ratchet, i.e. a periodic structure with broken spa- 
tial symmetry, yields the possibility to obtain a directed 
current in the presence of noise and unbiased nonequilib- 
rium forces (!]]. As such, ratchets belong to the class of 
so termed Brownian motors being devices operating 
far from equilibrium, and which combine noise and asym- 
metry to generate particle's transport. A huge amount 
of experimental and theoretical work exists on ratchets 
ruled by classical thermal fluctuations jl], |^| , also due to 
their possible application in biological systems ||, Q. In 
contrast, little is known about the ratchet effect in the 
quantum realm. This is partly due to the challenge in the 
experimental realization of suitable ratchet potentials. 
Only recently rectification of quantum fluctuations in tri- 
angularly shaped semiconductor heterostructures S , and 
in quasi one-dimensional Josephson junction arrays || 
has been observed. Theoretically, this originates from the 
complexity in the investigation of the quantum ratchet 
dynamics. By now the quantum ratchet effect has been 
tackled only for adiabatically rocked ratchets in a mostly 
numerical work fij, for peculiar single band models ||, 
in the presence of an external force of on-off type M , and 
for a weakly corrugated potential [flO| . Typical quantum 
features, as e.g. a current inversion upon temperature 
decrease were predicted and demonstrated ||. So far, 
the band structure of the potential was not considered. 

In this letter for the first time a microscopic theory 
for the interplay among tunneling, vibrational relaxation 
and non-adiabatic driving in ratchet potentials sustain- 
ing few bands below the barrier, cf. Fig. 1, is presented. 
For these potentials the semiclassical requirement |7) of 
having many bands below the barrier is not met. Our 
treatment, mostly analytical, is based on the real time 
path-integral formalism for open quantum systems [ pT| . 
Noticeably, in the temperature and driving regime in 
which the dynamics is effectively restricted to the low- 
est band of the periodic potential, no current rectifica- 
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FIG. 1: Ratchet potential with three bands below the barrier 
(shaded regions). The potential height is in units of the dis- 
tance Tiuio between the centers of the second and first band, 
and the length is in units of the period L. 



tion occurs. In fact, a reduction to the lowest band of 
the ratchet potential only retains information about the 
periodicity of the original Hamiltonian, but not about 
its reflection properties. To take into account the vibra- 
tional motion within the well, and hence the asymme- 
try of the ratchet potential leading to the ratchet effect, 
at least two bands should contribute. Starting from the 
system-plus-reservoir model of a driven ratchet system 
bi-linearly coupled to a harmonic oscillator bath, we pro- 
ceed in three steps: i) We consider the regime of tempera- 
ture and driving parameters which allows the truncation 
of the full dynamics to the Hilbert space spanned by the 
M lowest bands of the periodic ratchet potential, ii) A 
rotation to the so-called discrete variable representation 
(DVR) Q , being the eigenbasis of the discrete position 
operator coupling to the bath, is performed. In this ba- 
sis, the isolated ratchet problem is described by a tight 
binding model with a periodicity of M sites, and with 
further nearest neighours coupling, iii) The bath degrees 
of freedom are traced out exactly in the DVR basis, and 
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a closed form expression for the averaged current in the 
incoherent tunneling regime is obtained. Subsequently 
our results are applied to investigate current rectification 
and its optimisation in ac-driven ratchets. 

We start considering the total Hamiltonian H(t) = 
H R + H cxt (t) + H B , where H K = p 2 /2m + V K (x) is 
the Hamiltonian for a particle of mass m moving in 
the asymmetric periodic ratchet potential Vr(x + L) = 
Vr(x), cf. Fig. 1. Due to Bloch theorem, the spec- 
trum of the isolated system follows from the Schrodinger 
equation Hn\^ nt k) — £ n (k)\^ nt k), with n denoting the 
band index, and with k being the wave-vector. Time- 
reversal symmetry plus periodicity yield £ n (k) — E n + 
Em=i(^«'/2) cos(mfcL). The action of the unbiased 
force is captured by H ext (t) = xF(t), with F(t) = 
F cos Qt, while Hb is the standard Hamiltonian for a har- 
monic oscillator bath with bilinear coupling in the bath 
and system coordinates . The character of the bosonic 
bath is then fully described by the spectral function J(ui). 
Such a quantity is obtained from the equilibrium noise 
properties of the real system to be described. We assume 
the Ohmic form, relevant e.g. in quasi one-dimensional 
Josephson junction arrays, J{oS) = r)u>/{\ + {lo/u>b>) 2 ), 
corresponding to a white noise spectrum in the classical 
limit, with 7/ being the friction coefficent, and with lub> a 
bath cut-off frequency. 

We introduce Wannier states localized in cell j as 
\n,j) = f^ L dke- ik ^ L \^ n , k ). Due to lack 

of inversion symmetry, it holds (x\j, n) :— ^j^ n (x) ^ 
^ r _j l „(— x) for some x. In the following, we restrict 
our attention to the lowest M bands of the poten- 
tial. We focus on Hamiltonians whose spectrum is well 
approximated by tight-binding Hamiltonians with only 
nearest- neighbors coupling (miniband-approximation) , 
i.e., £ n (k) = E n + (A„/2)cosfcL. 

Then, the ratchet Hamiltonian is fully characterized by 
the 2M parameters, {A n ,E n }, n = 1, ..,M, determining 
the band widths, and the band centers, respectively. It 
has the form in the Wannier representation 

oo M . 

H * = E E^>w+M+ij+M>o>i) 

j — — oo n—1 

oo M 

+ E E^i-?»o>i ■ (!) 

j=—oo n=l 

Correspondingly, the discrete position operator reads 

oo M 

a;= E E Xj(n,m)\j,n){j,m\ , (2) 

j=— oo n,m=l 

where, due to the localized character of the Wannier 
states, intercell matrix elements were neglected. Here, 
Xj(n, m) = {j, m\x\j, n) := jLS nm + £ nm , where the ma- 
trix £„ m = £ mn carries information on the shape of the 



ratchet potential Vr(x) within a cell. Only within the 
commonly performed lowest band truncation, M = 1, 
the position operator x is diagonal in the Wannier basis. 
We wish to evaluate the average particle's velocity 

v = <w M (t))n := ( lim Ti{xp(t)}) n . (3) 

t — >oc 

Here pit) = TrBathW^) is the reduced density matrix 
(RDM) of the system. It is obtained by performing the 
trace over the bath modes in the density matrix W(t) 
of the system-plus-reservoir. Finally, (...}o denotes the 
time average over a driving period. The evaluation of 
the velocity v is quite intricate: As clearly seen from 
(||) , the coupling to the thermal bath and to the external 
field modifies the quantum coherent Bloch picture in (Q). 
Being the position operator not diagonal in the Wannier 
basis, intcrband transitions become possible due both to 
the coupling to the thermal bath and to the external field. 

For a harmonic bath bilinearly coupled to the system, 
the trace over the bath degrees of freedom can be done ex- 
actly by use of real-time path-integral techniques. This 
requires, however, to express the full system-plus-bath 
Hamiltonian H(t) in the basis of the eigenstates of the 
position operator, the so termed discrete variable rep- 
resentation (DVR) @ The contributions H cxt {t) 
and Hb to H (t) are already diagonal in this basis, but 
not the isolated ratchet Hamiltonian. Diagonalization of 
the matrix £ nm yields the orthogonal transformation U 
from the Wannier to the DVR representation, \j, n) = 
Yin ^n,/ib'iA*}j as wei l as the position operator eigenval- 
ues and eigenvectors x\j, p) = Xj 4L \j, p) = {jL + x^)\j, fi). 
Here /z labels the DVR states within a cell. Thus the 
ratchet Hamiltonian reads in the DVR basis 

oo / M M '' 

H *= E E £ mIj^)0^I+ E A^; a lj,M)(j,H 

J=-oc \fJ.=l u^v=\ f 

oo M 

+ E E A "7(|J,M>(.7 + M + \j + i,m)0>I) , 

j — — OO fl,U—l 

where the on-site energies read e u = Yn^riU n ^U niJl , 
while AJf^J a = J2 n EnUn^Unu are related to the miracell 
vibrational motion. Finally, A^ cr ee E„(A n /4)?7 njtl {7 nl/ 
accounts for intercell tunneling. Tunneling transitions to 
any DVR state in the neighboring cell are allowed. Ad- 
ditionally, because H cxt (t) is diagonal in the DVR basis, 
the Hamiltonian H(t) has time-dependent on-site ener- 
gies F(t)xj^ which add to the intrinsic ones e u . We have 
now all the ingredients to evaluate the asymptotic tun- 
neling current. Its expression reads 

oo M 

v as {t)=\im V y^Xj^Pj^t) , (4) 

t — >oo * — ' * — ' 

j=-oo u=l 

with the Pj,f_i{t) being the diagonal elements of the re- 
duced density matrix in the DVR basis. They repre- 
sent the probability of finding the particle in the state 
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\j, fi) at time t. In the following, we restrict our at- 
tention to the incoherent tunneling regime, reached at 
high enough temperatures and/or friction pj| . More- 
over, we consider moderate-to- high frequencies f2, such 
that uj > Q > r£j t , with T 3 ) 3 being the field-averaged 
transition rate from DVR site \j, v) to \j' , n). The upper 
bound, with Hluq — E2 — E\, is needed because of con- 
sistency with the truncation of the original problem to 
the first M bands. The lower bound enables to approxi- 
mate the evolution of the averaged occupation probabil- 
Pj,n(t))n by the coupled equations 



ities Pj, ^{t) : 



dt 



= -E[ r if^>(*)- r ^>(*)] 



v 

Ef^^-mC*) + nt^pj+iA*)] • ( 5 ) 



Within an incoherent tunneling description it is appro- 
priate to evaluate such rates as an expansion in terms of 
the coupling matrix elements A^, with A 3 ' 3 V — A^ lt t/ ra 
and Aj;^ 1 = A^ cr . To lowest order we find ' 

A;-! h) / dre-^W 

J —OO 

'* ^' *(t))^« 

where the external field parameters enter the argu- 
ment of the zero-order Bessel function Jo. This yields 
a strong reduction of the rate (|^), when — e„ = 
nhQ and the ratio F(xj tfl — xj/ „)/ftfi hits a zero of 
the Bessel function J n . In the absence of dissipa- 
tion, similar conditions would yield a dynamical local- 
ization of the particle 14, |l5}| . However, finite temper- 
atures of the bath hinder this effect. The environmen- 
tal influence is determined by the bath correlation func- 

tion Q]£(t) = {X] _\2'- r / °° rfc^[coth(ft/W2)(l - 
cos wt)-H sin wt] jn). Note that the dissipation strength 
depends on the square distance between the DVR states 
involved. It is position-dependent through the dimension- 
less coupling parameter 



2nh 



(7) 



For symmetric cosine shaped potentials a bath-induced 
localization at zero-temperatures has been predicted for 
dissipation parameters a = nL 2 /(2irh) > 1 p6| . A simi- 
lar tendence to localization occurs in our ratchet for large 
dissipation strengths ct 3 } 3 ^ 1 > 1 and zero temperature. 
However, appropriately tuned ac-fields can help prevent- 
ing localization. 

The stationary solutions of the coupled system of equa- 
tions can be found upon Laplace transformation. 



Then, the current assumes the compact form 



(8) 



where we introduced the forward and backward time- 

ri'j-F 1 , and the vibrational 



averaged rates T^v T ' 
relaxation rates T 



intra 



r 3 > 3 u (fx ^ v). The asymptotic 
cell occupation probabilities := lim^oo J^ . P JiM (£) 
are expressed in terms of the averaged rates. For M — 
1, the current in (^) is zero, as a consequence of time- 
reversal symmetry. Thus single-band quantum ratchets 
support no current, as recently reported in They 
define a new class of supersymmetric Q potentials in the 
quantum regime. At least two bands should contribute 
to the dynamics in order to have a non- vanishing current. 
For M = 3 contributing bands (or DVR states) we find 



P? 



r 2 r 3 



-^2,3-^3,2 



S/x^i/^y,^, and T u ^ 



(9) 



"pintcr,/ 1 pinter,b 



where T p 

r l J 1 ^ ra . Likewise, pSf and pf can be obtained from 
upon cyclic index permutation. In the absence of ac- 
driving the rates (|(3]) obey the detailed-balance relation 



F 3 ' 3 e~ 



(10) 



so that the total current in (||) adds up to zero. The 
ac-field invalidates (|To|), and a net current is originated. 

Up to this point, our results are general. We apply 
now our theory to make predictions for vortex motion 
in quantum ratchets devices. These could be realized 
with rectangular |(| or circular fl9[| arrays of different- 
sized Josephson junctions. Upon applying a magnetic 
field perpendicular to the array, current vortices are orig- 
inated, whose dynamics is homologous to that of a quan- 
tum Brownian particle in a periodic ratchet potential 
|p7[ @. The particle's mass is m — ^C/2a 2 , where $0 
is the flux quantum, C is the average capacitance of the 
junctions, and a is the average junction distance in the 
direction of motion. The magnitude of the potential felt 
by the vortices is proportional to the Josephson energy 
Ej. The vortices can be put into motion upon injecting 
a current / = Iac cos(f2t) into the array, which exerts a 
Lorentz-like force on such particles. At low vortex densi- 
ties, the voltage drop V across the array is related to the 
vortex velocity by the second Josephson relation 



V = sv$ /L 



(11) 



where s is the one-dimensional vortex density. Figures 2 
and 3 illustrate our predictions for current rectification 
for the driven ratchet potential shown in Fig. 1. It sup- 
ports only three bands below the barrier, being charac- 
terized by Ai = -4.43 x lQ- 3 hw , A 2 = 6.68 x I0~ 2 fuv , 
and A3 = — 3.10x 10 _1 ?iwo. The center of the bands read, 
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FIG. 2: Output velocity (or dc-voltage for the vortex ratchet 
case) vs. amplitude of the driving ac-signal (ac-current for 
vortex ratchet) for different driving frequencies Q. The di- 
mensionless friction parameter a = 2-KTjL 2 /h is choosen to be 
a = 0.5, describing moderate damping. Negative or positive 
ratchet currents are obtained upon appropriate tuning of the 
driving field parameters. 
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the current on the applied ac-amplitude and frequency. 
Likewise a non-monotonic behavior occurs varying damp- 
ing. For example, different current directions are ob- 
served for small ac-amplitudes in the curves with driving 
£1 = ujq of Fig. 2 and Fig. 3. This inversion also oc- 
curs in combination with a change in the bath cut-off 
frequency wd ■ As damping is further increased the main 
resonance peaks get shifted to even larger values of the 
driving amplitude, and the ratchet signal is reduced. 

In conclusion we investigated the ratchet mechanism 
in periodic quantum structures with few bands below 
the barrier, and driven by nora-adiabatic time-dependent 
fields. A set of coupled master equations for the re- 
duced density matrix elements in the basis which di- 
agonalizes the position operator are derived. The va- 
lidity of our method is restricted to the deep quantum 
regime, since thermally activated processes over the bar- 
rier are disregarded. Our approach is complementary to 
the semiclassical analysis of M , where many bands below 
the barrier are required. We discussed current rectifica- 
tion in Josephson ratchets. Current reversals result from 
the interplay of driving-induced dynamical-localization 
|p^| versus dissipation-assisted diffusion. Optimal cur- 
rent rectification is obtained for moderate damping, and 
driving frequencies of the order of the typical relaxation 
time of the environment. 

We thank M. Thorwart, P. Hanggi and J.E. Mooij for 
useful discussions. Support by FOM is acknowledged. 



FIG. 3: Same as in Fig. 2 for larger friction, a — 2.0. The 
main resonances occur at larger driving amplitudes. 



Ex = 5.87 x lO^ftwo, E 2 = 1.59fiw , and E 3 = 2A5hiu . 
This is the potential reported in [p] at a temperature 
of T — 50mK, corresponding to a charging energy 
Ej = 429fj,eV—8hujQ. We choose an average capacitance 
C = 2fF, an average distance a = L/2> = 2/im, and frus- 
tation s = 0.277. The average junction resistance R is 
taken as a free parameter. It determines the friction pa- 
rameter 7/ w 2 ^° 2 , and the bath cut-off uj^ « (i?C) _1 
p8| . In Figs. 2 and 3 we focus on two different val- 
ues of the dimensionless friction parameter a, cf. below 
(0), i.e., a = 0.5 and a = 2.0 which describe a mod- 
erate and a strong damping regime, respectively. These 
fully specifies the ratchet potential and the bath param- 
eters. Finally, the injected ac-current induces a force 
F(t) = Iac cos(fit)<I>o/iV c a, where iV c = 304 is the chosen 
number of junction columns of the array. Figs. 2 and 3 
reveal how the interplay among driving-induced localiza- 
tion near the zeros of the Bessel function in (^), driving- 
assisted tunneling when the driving frequency matches 
resonances of the dissipative ratchet system, and dissipa- 
tion results in a nontrivial non-monotonic dependence of 
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